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Abstract 

With high energy heavy ion colhsions one tries to create a new 
forms of matter that is similar to the one present at the birth of our 
Universe. Recent development on flow pattern, initial energy-density 
and freeze-out temperature shows that most likely this new form of 
matter is in a deconfined state, has colored degrees of freedom and is 
more fluid-like than gas-like. In present paper we calculate estimations 
on the physical properties of this new-old matter. 

"We simply do not yet know enough about the physics of elementary 
particles to be able to calculate the properties of such a melange with 
any confidence. . . . Thus our ignorance of microscopic physics stands 
as a veil, obscuring our view of the very beginning. " 

S. Weinberg, about the first hundredth of a second [1] 



1 Introduction 

Ultra-relativistic collisions, so called "Little Bangs" of almost fully ionized 
Au atoms are observed in four major experiments at the RHIC accelerator at 
the highest currently available colliding energies of ^/sNN = 200 GeV. The 
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aim of these experiments is to create new forms of matter that existed in 
Nature a few microseconds after the Big Bang, the creation of our Universe. 

Quantum Chromodynamics (QCD), the theory of quarks and gluons, 
the strong force interacting between them and their color degree of free- 
dom was formulated and established soon after Weinberg's famous book |1] 
about the early Universe had been published. Confinement is an important 
(though mathematically never proven) property of QCD, its consequence is 
that quarks are bound into hadrons in a matter of normal temperature and 
pressure. 

In the early Universe, energy density was many orders of magnitude 
higher than today, and at that high energy densities, deconfined phases of 
colored matter might have existed. Quark-gluon plasma (QGP) is such a 
phase, that might have existed during the first few microseconds after the 
Universe came into existence. This type of matter was searched for at the 
SPS, and experiments at RHIC are continuing this effort. Evidence for for- 
mation of a hot and dense medium in gold-gold collisions was found based 
on a phenomenon called jet quenching, and confirmed by its disappearance 
in deuteron-gold collisions [2]. 

A consistent picture emerged after the first three years of running the 
RHIC experiment: quarks and gluons indeed become deconfined, but also 
behave collectively, hence this hot matter acts like a liquid [3], not like an 
ideal gas theorists had anticipated when defining the term QGP. The situa- 
tion is similar to as if prisoners (quarks and gluons confined in hadrons) have 
broken out of their cells at nearly the same time, but they find themselves 
on the crowded jail-yard coupled with all the other escapees. This strong 
coupling is exactly what happens in a liquid [1]. 

1.1 A sign for hydrodynamic behavior: elliptic flow 

Azimuthal asymmetry of single particle spectra measured in relativistic heavy 
ion collisions is called elliptic flow (f2). It is an indication of liquid-like be- 
havior [S], and can be explained by hydrodynamics [HI El ISj- In the hydro- 
dynamic picture it turns out, that elliptic flow can result from the initial 
spatial asymmetry but also from momentum-space asymmetry. Important 
is, that in contrast to a uniform distribution of particles expected in a gas- 
like system, this liquid behavior means that the interaction in the medium 
of these copiously produced particles is rather strong, as one expects from 
a fluid. Detailed investigation of these phenomena suggests that this liquid 



2 



flows with almost no viscosity [H]. 



1.2 Relativistic perfect fluids 

Perfect hydrodynamics is based only on local conservation of charge and 
energy-momentum and on the assumption of local thermal equilibrium, and 
this is the tool that we use to describe and calculate the properties of the 
matter created in relativistic heavy ion collisions at RHIC. While there are 
accelerating non-relativistic solutions in the literature, recent development 
shows also relativistic solutions that can be compared to the data [lOl [TT] . 
Local conservation of charge and four-momentum reads as 

a,(nM") = 0, (1) 
d^T^"" = 0, (2) 
T^"' = (e + p)u^'u'' - pg^'r (3) 

We find the following solution (for arbitrary Xind=l, k=1 and for 
A = 2 in arbitrary d with k = d) [TU] : 



V 



tanh Xrj (4) 



Ad 



n = no [—j z/(s), (5) 
T = Top)"'" (6) 

where v{s) is an arbitrary function of the scale variable s, and rj is the pseudo- 
rapidity. As an illustration, fluid trajectories of this solution are shown on 
fig. [U See details in ref. |10j . 



2 Results 

2.1 An advanced estimate on the initial energy density 

Based on the above solution of eqs. 0116] let us estimate the initial energy den- 
sity of relativistic heavy ion or p-|-p reactions. As our solution is an acceler- 
ating one, and we do not neglect the initial acceleration period, we improve 
the renowned Bjorken estimate both quantitatively and qualitatively. Simi- 
larly to Bjorken's method [12], we can estimate the initial energy density (see 
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Figure 1: Fluid trajectories of the new exact solution of perfect fluid hydro, 
corresponding to d = 1 and A = 2. The trajectories are shown both inside 
and outside the lightcone. 

ref. [TU] for details). Finally we get a correction to the widely used Bjorken 
formula, depending on the acceleration parameter (A): 



The acceleration parameter can be extracted from the measured rapidity 
distribution [10]. For fiat rapidity distributions, a ^ oo (A — > 1, i.e. no ac- 
celeration) and the Bjorken estimate is recovered. For A > 1, the correction 
factor is bigger than 1. Hence we conclude that the initial energy densi- 
ties are under-estimated by the Bjorken formula. For realistic RHIC data 
from BRAHMS [13], the correction factor can be as big as s/sBj ~ 2.2 [TO] . 
Thus smaller initial bombarding (or colliding) energies are needed to reach 
the critical energy density in high energy heavy ion collisions, than thought 
previously using Bjorken's renowned formula. 

2.2 Estimating the freeze-out temperature 

We estimated the freeze-out temperature of these Little Bangs, fitting data 
to the Buda-Lund hydro model [HI [II]. Recently, Fodor and Katz calculated 
the phase diagram of lattice QCD at finite net baryon density: their results 
indicate that the transition from confined to deconfined matter is a cross-over 
with a nearly constant critical temperature, Tc = 175±2 MeV [15]. The result 
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Figure 2: Left: charged particle ^ distributions of ref. [13] fitted with the 
result of the relativistic hydro solution of ref. [10]. Right: the correction 
parameter (obtained from fits show on the left panel) as a function of freeze- 
out time versus thermalization time (tj/to). At reasonable values of 10-15, 
the correction to the Bjorken estimate of energy density is a factor of ~2. 



of the Buda-Lund fits to RHIC Au+Au data of refs. [131 [ISl [HI [Tg Hg [201 [2ll 
[22] (shown on fig. [3]), in particular the value of the fit parameter Tq (central 
freeze-out temperature, see details in ref. [23]), indicates the existence of 
a region several standard deviations hotter than the critical temperature. 
This is an indication on quark deconfinement in Au -|- Au collisions with 
= 130 and 200 GeV at RHIC [231 [2l [25], confirmed by the analysis 
of pt and 7] dependence of the elliptic flow [6] . A similar analysis of Pb+Pb 
collisions at CERN SPS energies yields central temperatures lower than the 
critical value, Tq < [221 [27] . 



2.3 Universal scaling of the elliptic flow 

The Buda-Lund calculation of the elliptic flow results (under certain condi- 
tions detailed in ref. [6]) in the following simple universal scaling law: 

^2 = fH, 8 
Io{w) 

thus the model predicts a universal scaling: every V2 measurement is pre- 
dicted to fall on the same universal scaling curve Ii/Iq when plotted against 
the scaling variable w (see details in ref. [28]). 
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BudaLund hydro fits to 130 AGeV Au+Au 



BudaLund hydro fits to 130 AGeV Au+Au 
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Figure 3: The upper four panels show a simultaneous Buda-Lund fit to 0- 
5(6) % central Au+Au data on pt and r] spectra and HBT radii at ^/suN = 130 
GeV [Ml [HI [m [191 ED] . The lower four panels show similar fits to 0-30 % 
central Au+Au data at a/snn = 200 GeV [131 1211 [22]. Fit parameters are 
summarized in ref. 1231. 
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This means, that V2 depends on any physical parameter (transverse or 
longitudinal momentum, center of mass energy, centrality, type of the collid- 
ing nucleus etc.) only through the (universal) scaling parameter w. 

In ref. [28] we have shown that the excitation function of the transverse 
momentum and pseudorapidity dependence of the elliptic flow in Au+Au 
collisions (RHIC data from refs. [SllSHlEn]) is well described with the formulas 
that are predicted by the Buda-Lund type of hydrodynamical calculations. 
We have provided a quantitative evidence of the validity of the perfect fluid 
picture of soft particle production in Au+Au collisions at RHIC up to 1- 
1.5 GeV but also show here that this perfect fluid extends far away from 
mid-rapidity, up to a pseudorapidity of r^beam — 0.5. The universal scaling of 
PHOBOS V2{'n) [21], PHENIX V2{pt) [S] and STAR V2{pt) [30], expressed by 
Eq. ([8]) and illustrated by Fig. Hie provides a successful quantitative as well 
as qualitative test for the appearance of a perfect fluid in Au+Au collisions 
at various colliding energies at RHIC. 

2.4 Chiral symmetry restoration 

Correlation functions are important to see the collective properties of parti- 
cles and the space-time structure of the emitting source, e.g. the observed size 
of a system can be measured by two-particle Bose- Einstein correlations [31] . 
The mt dependent strength of two-pion correlations, the so-called A,, param- 
eter, which is related to the extrapolated value of the correlation function at 
zero relative momentum, can be used to extract information on the mass- 
reduction of the ?7' meson, a signal of Ua{^) symmetry restoration in the 
source [321 [33l [Ml [35] . 

PHENIX analyzed [SB] A^.(mt) with fits to two-pion correlation func- 
tions using three different shapes. Gauss, Levy and Edgeworth, described 
in refs. [271 [361 [37] . A comparison of the measurements with model calcula- 
tions of ref. [32] using FRITIOF results for the composition of the long-lived 
resonances and a variation of the rj' mass is presented in fig. [5l If we re-norm 
the A*(mj) curves with their maximal value on the investigated rrit interval, 
they overlap, confirming the existence and characteristics of the hole in the 
A*(mt) distribution. 

Gauss fit results agree with former PHENIX measurements (see ref. [22]). 
Regarding Ua{^) symmetry restoration, conclusion is that at present, results 
are critically dependent on the understanding of statistical and systematic 
errors, and additional analysis is required to make a definitive statement. 



7 




Figure 4: PHOBOS [29] (a-b), PHENIX [5] (c) and STAR [30] (d) data on 
elliptic flow, 172, plotted versus pt and rj and fitted with Buda-Lund model. 
Elliptic flow versus variable w is plotted in panel (e): data points of plots 
(a-d) show the predicted [6] universal scaling. See fit parameters in ref. [28] 

3 Summary and conclusions 

In summary, we can make the definitive statement, based on elliptic flow 
measurements and the broad range success of analytic hydro models, that in 
relativistic Au+Au collisions observed at RHIC we see a perfect fluid. Based 
on our estimates on the temperature and energy density we also conclude 
that the observed matter is in a deconfined state. We also see a possible 
signal of partial symmetry restoration in the mass reduction of r/' bosons. 
Future plan is to explore all properties of the Quark Matter, by analyzing 



8 



X (If If ) - RUN4 200GeV Au+Au 



X (jf If) - RUN4 200GeV Au+Au 



X (It* It*) - RUN4 200GeV Au+Au 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




PHENIX PRELIMINARY 



= 600 MeV 
= 500 MeV 
= 400 MeV 
= 300 MeV 
= 200 MeV 




i Gauss fit 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 p 
1.1 



1 
0.8 
0.6 
0.4 
0.2 





PHENIX PRELIMINARY 




Levy fit 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9, 



0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




PHENIX PRELIMINARY 



- m.„. = 958 MeV 
m.;. = 600 MeV 
m. , = 500 MeV 
m.;, = 400 MeV 

- m. . = 300 MeV 

- m.;, = 200 MeV 



♦ Edgewortli fit 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9, 



1.0 



0.9 



0.8 



0.7 



PHENIX PRE 



LIMINARY 



Run4 200GeV Au+Au 



3 Edgeworth (k^=0.1 fixed) 

H Levy (a = 1.35 fixed) 
<► Gauss 



0.6 

0.20 0.25 0.30 0.35 ^ 0.40 0.45 0.50 0.55 

t 

Figure 5: Top three figures: measured K{mt) of ref. [36] compared to calcula- 
tions using the model of ref. [32] with various rj' mass values. Bottom figure: 
A*(mt) curves with equal number of fit parameters and re-normed with their 
maximal value on the interval of 0.20GeV< rrit <0.55GeV all show the same 
shape. See details of the fits in ref. [36] . 



more data and using higher luminosity. We are after the full map of the 
QCD phase diagram, and in order to explore it, we also have to go to higher 
energies and compare them to lower energy data. If the Quark Matter is the 
New World, then Columbus just realized he is not in India, but on a new 
continent. 
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"It does not make any difference how beautiful your guess is. It does 
not make any difference how smart you are, who made the guess, or 
what his name is — if it disagrees with experiment it is wrong. " 

R. P. Feynman, about discovering new laws [38] 
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